SO (3) symmetry between Neel and ferromagnetic order parameters for graphene in a 
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I consider the Hubbard model of graphene in an external magnetic field and in the Hartree-Fock 
approximation. In the continuum limit, the ground state energy at half filling becomes nearly 
symmetric under rotations of the three-component vector A = (Ni,N2,m), with the first two 
components representing the Neel order parameter orthogonal to and the third component the 
magnetization parallel with the external magnetic field. When the symmetry breaking effects arising 
from the lattice, Zeeman coupling, and higher Landau levels are included the system develops a 
quantum critical point at which the antiferromagnetic order disappears and the magnetization has 
a kink. The observed incompressible states at filling factors ±1 are argued to arise due to a finite 
third component of the Neel order parameter at these electron densities. Recent experiments appear 
consistent with Ni = N2 = 0, and N3 7^ 0, at the filling factors zero and one, respectively. 



I. INTRODUCTION 



The nature of the ground state of two-dimensional car- 
bon (i. e. "graphene") in a uniform magnetic field is 
an issue which has attracted a lot of attention lately. 
The interest mainly stems from the recent experimen- 
tal observation [l[ of the additional plateaus in the Hall 
conductivity at the values of zero and unity which are 
not naturally explained by the picture of non-interacting 
electrons alone Q, Q, Q. The mechanism behind the 
formation of the gaps in the energy spectrum that appar- 
ently develop within the zeroth Landau level has been a 
matter of debate. On one hand, the Coulomb repulsion 
may on generalgrounds be expected to favor breaking of 
the sublattice (5J , or "valley" symmetry Q , Q , [H , which 
in the zeroth Landau level are equivalent. If the Zeeman 
coupling of the electron spin to the magnetic field is en- 
tirely neglected the sense of this symmetry breaking can 
be equal or opposite for the two projections of spin, de- 
pending on details of the interaction on the lattice scale. 
This would yield a staggered pattern of either charge or 
spin in the ground state, with the concomitant many- 
body gap in the quasiparticle spectrum Q. One could, 
however, also imagine the opposite limit of a strong cou- 
pling of the electron spin to the magnetic field, in which 
the gap in the spectrum at half filling would become es- 
sentially a single-particle Zeeman gap. Of course, such 
a ferromagnetic ground state is already favored by the 
repulsive Coulomb interaction alone, due to the familiar 
physics of Hund's rule [lfj. The interplay between dif- 
ferent possible instabilities in graphene in magnetic field 
represents an important unsolved problem at the moment 
[Til ] . Here I would like to point out and explore a 

surprising symmetry between different order parameters 
which may help its resolution. 

I consider the minimal Hubbard model for the inter- 
acting electrons living on a honeycomb lattice. Although 
admittedly a simplification, the Hubbard model already 
contains most of the relevant physics. Without a mag- 
netic field it exhibits both the ferromagnetic and anti- 
ferromagnetic ground states in its phase diagram [l3T ]. 




FIG. 1: The dependence of the Neel order parameter (N±) 
orthogonal to and the magnetization (m.3) parallel with the 
external magnetic field on the Zeeman energy (A) at the fill- 
ing factor zero (Eqs. (15)-(18)). The critical point is located 
at A c = 2(g a — <?/)/(7r£ 2 ), where g a and <?/ are the coupling 
constants in the Neel and the ferromagnetic channels, respec- 
tively. I is the magnetic length. 



with the latter as stable at half filling [14J. In a field, 
constraining the Neel vector to lie parallel to it leads to 
a discontinuous transition between the antiferromagnetic 
and ferromagnetic ground states at a critical Zeeman cou- 
pling, as discussed in the previous work by the author . 
Here I calculate the ground state energy in the Hartree- 
Fock approximation, but for an arbitrary direction of the 
Neel vector. An immediate, and maybe not entirely un- 
expected result is that the ground state energy at half fill- 
ing is minimized when the Neel order parameter N is in 
fact orthogonal to the external magnetic field. Restrict- 
ing therefore N to its easy plane (N3 = 0) and taking the 
magnetization m to be parallel to the external magnetic 
field (mi = — 0), the ground state energy at half fill- 
ing becomes nearly symmetric with respect to an internal 
50(3) group of rotations of the three-dimensional order 
parameter A = (Ni, N2, m^). This symmetry becomes 
exact in the zeroth Landau level approximation, if the 
Zeeman effect and the discreteness of the lattice are also 
neglected. With these realistic symmetry breaking ef- 
fects included one finds that with the increase of Zeeman 
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coupling the system at half filling suffers a continuous 
quantum phase transition from the "mixed" state with 
both N± and 7713 ^ to a pure ferromagnet with 
N± = and m 3 7^ 0. This situation is depicted in Fig. f . 
With a change in the chemical potential the third Neel 
component N 3 eventually becomes finite and produces 
additional incompressible states at the filling factors ±1. 
An observation of such, albeit according to our estimate, 
a rather weak Neel order along the magnetic field at ±1 
would provide a direct support for our theory. 

Motivated by a recent experiment [l5j, I examine the 
dependence of the activation gap at filling ±1 on the in- 
plane component of the magnetic field, at the fixed per- 
pendicular field. The present theory predicts such a gap 
to be completely independent of the in-plane component 
only if at the filling zero N±_ = and 1113 7^ 0, in which 
case it is simply given by 2N 3 . Given that this seems to 
be the situation in the experiment, we conclude that for 
the experimentally relevant values of the parameters the 
ground state should have a finite A3 at the filling fac- 
tor ±1, whereas it should be a pure ferromagnet at zero 
filling. This conclusion is in accord with the observed 
dissipative nature of the state at the filling zero, which 
can be understood in terms of the particular edge-state 
transport implied by such a ferromagnetic ground state 
[l6| . It is then also in agreement with one of the possi- 
ble scenarios discussed previously in ref. (9), where the 
Neel order parameter was constrained to lie only along 
the magnetic field. 

The paper is organized as follows. In the next section 
I derive the expression for the ground-state energy of 
the system as a function of Neel and uniform magnetiza- 
tions, in the Hartree-Fock approximation. In the section 
III the emergence and the breaking of the approximate 
SO (3) symmetry by different terms in the Hamiltonian 
is discussed. The discussion of the result as they pertain 
to experiment is given in the section V. Finally, the main 
results are summarized in section VI. 



II. HARTREE-FOCK GROUND-STATE 
ENERGY 

Let me define the Hubbard model on a honeycomb 
lattice as H t + Hjj, where 

H t = -t J2 uUA)Va(A+bi)+H.c., (1) 



and 



H " = TSEIW^) + n ( A +b)) 2 + HA) - n(A + b)) 2 (2) 
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-(m(A) + m(A + b)) 2 - (m(A) - m(A + b)) 2 



sublattice is then at B = A + b, with b= (1/V3, l)(a/2). 
a is the lattice spacing, and n{A) = u\{A)u rT {A) and 
rh(A) = ul r (A)a„cr'U a '(A) are the particle number and 
the magnetization vector at site A. Variables at the sec- 
ond sublattice are analogously defined in terms of fermion 
operators v a (B). It is easy to check that Hjj is just the 
standard Hubbard on-site repulsion, written here in the 
rotationally invariant form which will prove to be suitable 
for our purposes. 

Here we are after the Hartree-Fock ground state of 
the Hubbard model in the uniform magnetic field and 
at half filling. Even at the highest laboratory fields 
~ 40T the magnetic length is much larger than the lat- 
tice spacing, I ^> a. It suffices therefore to consider 
only the continuum theory corresponding to the Hub- 
bard model. General form of the low-energy field the- 
ory of graphene has been discussed in detail before [Til ]. 
For present purposes it will be enough to consider the 
Hartree-Fock ground state obtained after the usual de- 
coupling of the interaction term in the ferromagnetic and 
antiferromagnetic channels. We thus assume a uniform 
density ((n(X)) = 1, X = A,B), and allow for both the 
uniform (m = (fh(A) + m(A + b)}) and staggered mag- 
netizations (N — (rh(A) — m(A + b))). Standard manip- 
ulations [171 ] give the ground state energy per unit area 
and at half-filling to be 



N 2 m 2 
E= — + — + E {N,m], 
4ga 4 5/ 



(3) 



The sites A denote one triangular sublattice of the honey- 
comb lattice, generated by linear combinations of the ba- 
sis vectors a\ = (V3, — l)(a/2), 0,2 = (0, a). The second 



where E§ is the ground state energy per unit area of the 
resulting single-particle Hamiltonian 

H H f = h <S> Hq - (N ■ a) <g 70 + ((A + m) • a) <g> I A . (4) 

X = gzB represents the Zeeman effect of the uniform 
magnetic field B, and gz is the electron g-f actor. Ho = 
i"/o"/i(~idi — Ai), with B = e^diAj, is the standard 
Dirac Hamiltonian near the two Fermi points in the spec- 
trum of H t [H[ , . For simplicity, the magnetic field 
will be assumed to be orthogonal to the graphene plane, 
until further notice. We work in "graphene representa- 
tion" of the Clifford algebra [H] in which 70 = I2 <£> 03 , 
7i = 03 <£> 02, 72 = —h ® 01, where {I2, <?} is the stan- 
dard Pauli basis of two-dimensional matrices. Likewise, 
J4 represents the four-dimensional unit matrix. In our 
units U = e/c = Vf = 1. The coupling constants gj 
and g a are both proportional to the original repulsion 
energy U , and positive. Although they are exactly equal 
at the lattice scale, as evident from the form of Hu in Eq. 
(2) , they in general will not be in the effective linearized 
Dirac theory sensible only below a certain momentum 
cutoff A. Subscribing to the usual logic of the low-energy 
description, <?/ and g a must be considered only as effec- 
tive parameters which themselves depend on the cutoff A 
in a way that ensures the cutoff independence of all phys- 
ical quantities. This observation will play a role in the 
selection of the ground state, as will be discussed shortly. 
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To proceed one needs the spectrum of the eight- 
dimensional Hamiltonian H^p. Perhaps this is most 
easily computed by casting it into a block-diagonal form 
H HF =H 1 ®H 2 (20|, with 

H 1(2 ) = ±h ® <r\{-idi - A\) -h® cr 2 (~id 2 ~ A 2 ) (5) 
-(A • <r) <8> ct 3 + (A + TO 3 )cr 3 h. 

Here we have chosen the third spin axis and the magneti- 
zation to be along the magnetic field, i. e. mi = mi = 0, 
rn 3 > 0. Next we note that Hi = UlHUx, with 
Ux = I2 @ i&2 , and 

H = ffo + «A r i7o73 + iA r 27o75-N37o + i(A + TO3)7375 ( 6 ) 

and 73 = 01 ® as, 75 = 02® 02. Similarly, H 2 = UlHU 2 , 
with U 2 — ia 2 Q)I 2 , with A3 — ► — A3. It is sufficient there- 
fore to find the spectrum of the single four-dimensional 
Hamiltonian H. After some algebra the eigenvalues of 
H up are this way found at ±e n(T with 

e„ CT = {A 2 + [(JVf + 2ns) 1 / 2 + <r(A + m 3 )] 2 } 1 / 2 , (7) 

with a = ±1, Aj_ = (Ni,N 2 ), and the degeneracies per 
unit area of 1/(tt1 2 ) and l/(27r/ 2 ), for n — 1,2,3, ... and 
n = 0, respectively. Note that for TV = the eigenvalues 
become the familiar relativistic Landau levels split by the 
Zeeman term, as expected. For A + rn 3 = 0, on the other 
hand, the spectrum reduces to the Landau levels of the 
massive Dirac Hamiltonian, with the mass |A|, also in 
agreement with the previous calculations [9j. 

At half filling then the ground state of Hjjf has all 
the eigenstates with negative energies filled, and all the 
others empty, so 



E [N,m 3 ] = 



1 



2tt/ 2 



(8) 



<T=±X 



The variational Hartree-Fock ground-state energy is then 
determined by the least value of the expression in Eq. (3). 



III. SO(3) SYMMETRY AND ITS BREAKING 

Let us minimize the Hartree-Fock ground state energy 
given by Eqs. (3) and (8), assuming A 7^ 0. Choosing \N\ 
and A3 as independent variables and then differentiating 
with respect to A3 yields 



<j(X + m 3 )[— + J2 



A3 



a=±X 



2 + 2nB 



0. (9) 



The left-hand side of the equation vanishes for A3 = 
0, and otherwise is a negative definite function of -/V 3 . 
7V3 = is therefore the only solution. One can also show 
that this solution represents the minimum of the energy. 
Restricting the Neel vector then to be orthogonal to the 



magnetic field the ground state energy in Eq. (3) can be 
rewritten as, 



E — ^so(3) + E i 



(10) 



where 



E 



Nl + m 2 1 
SO(3) -^^-^ 2 ' 



(Nl+ml)^, (11) 



and E' = E A + E z + E H ll, with 

p _ (ga-g/) „2 
h/\ = — : m 



3' 



(12) 



] ^ = ^2 + m D 1/2 ~ (#1 + ( A + ™ 3 ) 2 ) 1/2 ], (13) 



Ehll = ~ 2^2^ ^2 • 

er=±l n=l 

This form makes it manifest that the Hartree-Fock 
ground state energy at half filling is nearly symmetric 
with respect to rotations of the three-component order 
parameter A = (Ai, N 2 ,m 3 ). The identification of this 
approximate internal SO(3) symmetry is our central re- 
sult. The 5*0(3) symmetry is in the Hubbard model bro- 
ken by three terms of different physical origin: 1) the 
difference between the coupling constants g a and g m at 
the cutoff A (Ea), 2) the finite Zeeman coupling to the 
magnetic field A (Ez), and 3) the n ^ Landau level 
contribution (Ehll)- If one were to consider only the 
zeroth Landau level, use the bare values of the couplings 
for which g a = g mi and neglect the Zeeman coupling, 
the Hartree-Fock ground-state energy would become per- 
fectly SO(3) symmetric. In particular, the antiferromag- 
netic state with the Neel vector orthogonal to the exter- 
nal magnetic field and the ferromagnetic state with the 
magnetization along the same field would in this approx- 
imation appear as degenerate. 

Clearly, the Zeeman term favors magnetization. Simi- 
larly, it can be shown that Hull prefers the Neel compo- 
nents [2l|. It is less obvious what the sign of the "easy- 
axis anisotropy" , i. e. of g a — 5/, in Ea should be. This 
is determined by the flow of the two coupling constants 
as the high-energy modes between the momenta ~ 1/Z 
and ~ 1/a are integrated out essentially at zero mag- 
netic field. Since the leading instability of the Hubbard 
model on a honeycomb lattice and at half filling as the 
interaction strength is increased is towards the antiferro- 
magnetism, we may assume that in general g a > <?/ in 
the low energy theory. This would also be in accord with 
the explicit renormalization group calculation for weak 
couplings [14| . 

The ground state in the Hartree-Fock approximation 
is thus the result of the competition between the high- 
energy modes represented by H\ and Hull, and the 



4 



Zeeman coupling, which favor antiferromagnetic and fer- 
romagnetic solutions, respectively. Since the effect of 
Hhll is in the same direction as of H\, to keep the 
algebra simple we will drop Hhll altogether and just 
assume g a ~ 9f > in the Eq. (10). This is a partic- 
ularly good approximation at the laboratory magnetic 
fields at which the inclusion of the higher Landau levels 
is expected to provide corrections to our results of higher 
order in the small parameter a/l. Minimizing the energy 
in Eq. (10), for A < A c we then find 



\Nx\ = 



9a 



9a- gf 



91 



and for A > A r 



9a - 9f 



TV i =0, 



-A, 



m 



91 



9a~ 91 



(15) 



(16) 



(17) 



(18) 



where A c = 2(g a — gf)/(jrl 2 ). With the increase of the 
Zeeman coupling there is a continuous transition at which 
the Neel order disappears, and the magnetization satu- 
rates to its maximum (Fig. 1). 

Note that if the Neel vector is constrained to be along 
the magnetic field, so that N±_ = 0, the above quantum 
critical point becomes replaced with the level crossing be- 
tween the purely antiferromagnetic (iV 3 7^ 0, m 3 — 0) and 
purely ferromagnetic (N 3 = 0, m 3 7^ 0) ground states @. 
Such a pure antiferromagnet would correspond to a local 
maximum of the energy in Eq. (3), however. It seems al- 
ways energetically favorable to compromise between the 
two competing magnetic orderings by orienting the Neel 
vector orthogonally to the field. 



IV. DISCUSSION 

While it is not a priori clear at which side of the transi- 
tion should the experimental samples lie at half filling, in- 
creasing sufficiently the component of the magnetic field 
parallel to the graphene layer should always place the 
system into the purely ferromagnetic state. The exper- 
imental fact that the system at the filling factor zero is 
dissipative suggests that this is probably already the case 
even for the magnetic field completely orthogonal to the 
graphene's plane [liij]. Spin splitting as the origin of the 
measured gap at half filling has also been suggested very 
recently in ref. [l5|. Most importantly, the experiment 
of Jiang et al. offers compelling evidence that the gap 
at filling factor ±1 is independent of the total magnetic 
field, and thus likely to be a consequence of interactions 
[l5| . And this is precisely what follows from our theory 



@ . When the chemical potential gets close to the energy 
of the first excited states laying at 



ti c = ±y/Nl + (\ + m 3 )*, 



(19) 



the Hartree-Fock ground state energy may always be low- 
ered by developing a finite third Neel component N 3 . 
This pushes half of the states in question below the chem- 
ical potential and opens a gap. This way an incompress- 
ible state at filling factors ±1 would be formed irrespec- 
tively of whether N± was zero or finite at the filling factor 
zero. Note, however, that the Eq. (7) implies that this 
mechanism is operative only in the n — Landau level; 
the states belonging to the other Landau levels do not 
get split but only shifted in energy when N3 ^ 0. Con- 
sequently, at weak coupling iV 3 will become finite only 
in the ±1 state and otherwise not. The incompressible 
states that follow from the inclusion of a finite chemical 
potential into Eq. (3) lie therefore only at filling factors 
0, ±1, and all even integers. 

Let us assume therefore a week third Neel component 
iV 3 at the filling factor ±1. The activation gap at this 
filling then becomes 



K — 

^gap — 



2(A + m 3 ) 



(Nl + (A + m 3 ) 2 )i/2 



N 3 + 0(N$) 



(20) 



For a small Zeeman term, A < A c , Eqs. (15) and (16) 
give then 



E„ 



2A 



N 3 + 0(N 2 ), 



whereas when A > A r and Ni = 0. 



E, 



gap 



2N 3 + 0(Ni) 



(21) 



(22) 



Recalling that both N 3 ~ B± and A c ~ B±, whereas 
A ~ B, where B is the total and B± the perpendicular 
component of the magnetic field, we see that the gap is 
independent of the field's in-plane component only in the 
latter case. The experiment is therefore consistent with 
A^ 3 7^ at the filling factor ±1, and with N± = at the 
filling factor zero. This may not be very surprising in 
view of the equality between the couplings g a and gf at 
the lattice scale, which makes the critical Zeeman cou- 
pling Ar likely to be exceeded by the experimental value 
of A [H. 

Finally, restoring a finite range to electron-electron 
interactions may replace the antiferromagnetic ground 
state discussed here with a charge density wave with the 
particle density alternating (around the value of unity) 
on the two sublattices. In this case the mean field phase 
diagram would remain the same as discussed previously 
in ref. (9). Provided that the on-site repulsion is in- 
deed the dominant effect of the Coulomb interaction in 
graphene, our prediction is that at filing factors ±1 and 
at magnetic fields B ~ 10T the system has a weak Neel 
component N 3 ~ (a/l) 2 [iB ~ (10~ 3 — 10 _4 )/is per elec- 
tron. fiB is the Bohr magneton. The inherent weakness 
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of the Neel order derives from a small fraction of elec- 
trons occupying the relevant n = Landau level, which 
is the sole source of antiferromagnetism in this case. The 
existence of such ordering at the fillings ±1 can be used 
to distinguish the present proposal from all others in the 
current literature on the subject. 

V. SUMMARY 

To summarize, the ground-state energy in the Hartree- 
Fock approximation to the Hubbard model on honey- 
comb lattice and at half filling is nearly symmetric with 
respect to rotations between the components of the Neel 
order orthogonal to and the magnetization parallel with 
the external magnetic field. The effects of the symmetry 
breaking terms originating from discreetness of the lat- 
tice, Zeeman interaction, and n =/= Landau levels were 



examined. It was shown that whereas the component of 
the Neel order parameter parallel to the magnetic field 
always vanishes at half filling, it becomes finite (only) at 
filling factors ±1, introducing this way a gap in the spec- 
trum. The existence of such a weak Neel order along the 
magnetic field ~ 10 hb per electron is proposed as a 
litmus test of the present theory. 
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